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Recently synthesized 3D materials with Dirac spectrum exhibit peculiar electric transport qual¬ 
itatively different from its 2D analogue, graphene. Neglecting impuritiy scattering, the real part 
of the conductivity is strongly frequency dependent (linear), while the imaginary part is non-zero 
(unlike in undoped, clean graphene). The Coulomb interaction between electrons is unscreened as 
in a dielectric and hence is long range. We demonstrate that the interaction correction renders 
the electrodynamics nonlocal on a mesoscopic scale. The longitudinal conductivity <jl (related by 
charge conservation to the electric susceptibility) and the transverse conductivity <jt are different 
in the long wave length limit and consequently the standard local Ohm’s law description does not 
apply. This leads to several remarkable effects in transport and optical response. We predict a 
charging effect in DC transport that is a direct signature of the nonlocality. The optical response of 
the WSM is also sensitive to the nonlocality. In these materials p-polarized light generates bulk plas- 
mons as well as the transversal waves. The propagation inside the WSM is only slightly attenuated. 
At a specific (material parameter dependent) frequency the two modes coincide, a phenomenon 
impossible in a local medium. Remarkably, for any frequency there is an incident angle where total 
absorption occurs, turning the WSM opaque. 

PACS numbers: 

One of the common assumptions of electrodynamics 
in electrically active media is that the effect of exter¬ 
nal electric fields can be described locally by constitutive 
relations connecting the ” induced” currents to the elec¬ 
tric field even when spatial dispersion is present. Gener¬ 
ally, due to space - time translational symmetry of the 
material, the relation between Fourier components (cd 
is the frequency, k the wavevector) of the electric field 
and these of the induced current density within linear 
response reads: 

Ji (c j, k) = cTij (cu, k) Ej (u, k). (1) 

Here a is the AC conductivity tensor with indices i,j = 
x,y,z. The locality of the electrodynamic response in 
Fourier space means that the long wavelength limit ex¬ 
ists: dij (o;,k = 0 ) = (Jij (uj). The conductivity tensor 
for a homogeneous, isotropic, space and time-reversal in¬ 
variant (nongyrotropic) material at k = 0 simplifies into 
the simple form of Ohm’s law: 

(Tij (u) = Sijcr(uj) , ( 2 ) 

or J (cj) =g (lj) E (cj). 

On the microscopic level the locality is not 
guaranteed [lj. It hinges on the nature of the charges in 
the condensed matter system and presence of long-range 
interaction between them. These in turn determine the 
long-wave excitations of the material. In an insulator (or 
semiconductor at low temperatures) locality is simply a 
result of absence of gapless charged excitations. This 
does not apply to metals. 


In the free electron gas model of an ideal metal, i.e. 
neglecting both disorder and electron-electron interac¬ 
tions, there is no energy gap, so there are gapless charged 
excitations. The conductivity tensor in a metal can be 
uniquely decomposed into a transversal and a longitu¬ 
dinal part (assuming rotational and reflection symmetry 
for simplicity) [l|: 

^ a T (w, k) + ^-(t l (w, k), 

(3) 

Yet, the standard ’Lindhard’ type calculations in an ab¬ 
solutely clean metal with an arbitrary dispersion relation 
of the charge carriers and finite area Fermi surface, see 
Fig.la, show that the two scalars gt and gl are not in¬ 
dependent at small wavevectors, 

g t (cj,k = 0 ) = g l (u;,k = 0 ) = g (uj) ; (4) 
gl (cj, k) gt (td) k) = p (cd) k 2 + O (k A ) , 

thus leading to the local Ohm’s law Eq.(j2j). 

Impurities in a metal also define a length scale, the 
mean free path, that effectively makes the carriers’ mo¬ 
tion diffusive and the charge excitations ’’massive” in 
terms of their dispersion relation. Accounting for im¬ 
purities within the self consistent harmonic approxi¬ 
mation (yet still neglecting Coulomb interactions) thus 
makes the long wavelength limit of the Lindhard diagram 
smooth Beyond the harmonic approximation, com¬ 
posite excitations like diffusons have zero modes. How¬ 
ever coupling of the external fields to these excitations is 
’’soft” enough to cause the so-called ’’infrared divergen¬ 
cies” that in principle could make the long wavelength 
limit singular 


(cj,k) = 


ki kj 
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Returning to a very clean electron gas, the nonlocality 
can in principle arise due to long-range Coulomb interac¬ 
tions. However, long-range interactions are not screened 
only in insulators. In a metallic state with finite den¬ 
sity of charges (finite density of states on the Fermi sur¬ 
face) the nonlocality is prevented by the screening of the 
Coulomb force that becomes effectively short-range and 
thus unable to cause infrared divergencies. As a result 
clean interacting charged electron gas is local. 

All the above reasons ensuring locality, namely a di¬ 
rect energy gap, significant disorder and screening of the 
Coulomb interactions, are inapplicable to the recently 
discovered Weyl semi-metals (WSM). In these crystals 
electronic states are described by the Bloch wave func¬ 
tions, obeying the 3D ” pseudo-relativistic” Weyl (Dirac) 
equation with Fermi velocity v replacing the velocity of 
light, see Fig.lb. The novelty of the physics of WSM at 
Dirac point is clearly due to its ” ultrarelativistic” disper¬ 
sion, 5k = vk, of the elementary excitations. The Fermi 
level at neutrality point is occupied by a finite number 
of states (zero density of states), so that the excitations 
are those of the neutral plasma. 

Although the two-band electronic structure of bismuth 
was described by a four-component nearly massless Dirac 
fermion in 3D caused by spin-orbit interaction long ago 
[H (with spin replacing pseudospin), only recently several 
systems were experimentally demonstrated to exhibit the 
3D Dirac quasiparticles. Materials are quite diverse and 
include the time - reversal invariant topological Dirac 



), a bulk crys- 
CdsAs 2 with a 


semimetal ii Na^Bi (predicted 
tal symmetry protected semimetal 
single pair of Dirac points 0 , crystals on the phase 
transition boundary between topological and band insu¬ 
lators HgCdTe\ 11]. Recently, Weyl semi-metals exhibit¬ 
ing Fermi arcs on their surface were discovered 12-14). 
All these materials exhibit a great variety of new elec¬ 
tromagnetic transport and optical phenomena (not seen 
even in 2 D WSM like graphene) including giant diamag¬ 
netism, quantum magnetoresistance showing linear field 
dependence [HI Hsj . superconductivity [19] etc. 

There is no finite energy gap in WSM to ensure lo¬ 
cality as in band insulators or semiconductors. In this 
paper we assume that the WSM is clean enough, so that 
the disorder scale is irrelevant and the chemical poten¬ 
tial is tuned to the Dirac point. The elementary exci¬ 
tations are still ’’massless” at long wavelengths like in a 
metal. To no surprise the free (noninteracting) electron 
gas AC conductivity of WSM, calculated recently j 2 ()l.[ 2 lj|. 
turned out to remain local (see discussion below). How¬ 
ever calculating the screening due to the interaction 
corrections[ 20 ?, 22 ,^[, it became clear that the Coulomb 
interactions are still long range (namely like an insulator 
WSM is unable to screen the Coulomb interactions, un¬ 
like in the metal). The reason for locality mentioned with 
respect to metals therefore does not apply. So that there 
is no good reason to exclude the possibility that the long¬ 


wave electric response of WSM in the presence of strong 
interactions is nonlocal. To summarize WSM exhibits 
a curious mixture of properties usually attributed either 
to metals or insulators that does not allow to apply a 
conventional reasoning to establish locality of the linear 
response. 

In the present paper we demonstrate by explicit mi¬ 
croscopic calculation of the interaction effects in a rather 
generic model of WSM that the electrodynamics is indeed 
nonlocal. The rather unusual macroscopic electrodynam¬ 
ics then is formulated and applied to various physical 
phenomena. Several experimental setups in which the 
nonlocality can be demonstrated are suggested. These 
include the charging effects in DC transport, while an¬ 
other is the optical generation of both longitudinal (plas- 
mon) and transversal waves in these materials and their 
subsequent propagation. 


RESULTS 
The model 

An analogous calculation in graphene [26|,[27j (a 2 D ver¬ 
sion of WSM) reveals that in order to avoid complica¬ 
tions linked to the absence of scale separation in the sys¬ 
tem of relativistic massless fermions (known otherwise as 
’’anomaly”), one should use a well defined lattice model. 
Electrons in WSM are described sufficiently accurately 
for our purposes by the tight binding model of nearest 
neighbors on a cubic lattice [20u28j n = (see Fig.2a). 
The Hamiltonian is 

i+ hC ’ ( 5 ) 

n,i 

where cq are Pauli matrices, operators , a — 1,2, 
create a two - component spinor (describing two orbitals 
per site) and 7 is the hopping energy determining the 
Fermi velocity v = 'ya/h (a - lattice spacing). This de¬ 
fines a propagator depicted by arrows (see Fig.2b) in the 
Feynman diagrams. The Coulomb interaction (neglect¬ 
ing retardation effects) on the lattice can be viewed as 
an exchange of static photons. In term of Fourier com¬ 
ponents it takes the form (see[29] for details) 

v = ^ EE v p c kf p c kcflp<, (6) 

pkl 

where V is the sample volume and 

__ 

P sin 2 (p x a/ 2 ) + sin 2 (p y a/ 2 ) + sin 2 ( p z a/2 ) 

The photon propagator v p is represented by wavy lines 
in Fig.2. 
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It is convenient to calculate the transverse and lon¬ 
gitudinal conductivities from the scalar quantities using 
a l = kikj<Jij/k 2 and an = 2a t + 07 ,. The relations 
follow from Eq.(j3]). The calculation involves evaluation 
of diagrams 2c for free electrons and 2d, 2e, 2f for the 
leading interaction corrections, self energy, the vertex 
renormalization and ” glasses” diagram, respectively. De¬ 
tails involving renormalization, cancellation of both in¬ 
frared and ultraviolet divergences appear in the Supple¬ 


mentary Information SI1, while results are given Eq.(j8]) 
and Ea. (fTT|) . 

The electrodynamics of a free Fermi gas is local 

The dispersive conductivity tensor for a free (ne¬ 
glecting the electron-electron interaction) clean ultra- 
relativistic fermion gas at Dirac point is (see 0): 


(Jij (w, k) 


24 ttHvuj \ it s uj 2 


N is the number of Weyl fermions and A ~ 1/a is the ul¬ 
traviolet cutoff. In the long wavelength limit one recovers 
the AC conductivity of- (cj, k = 0) = ao (uj) Sij with 


CT 0 (uj) = 


Ne 2 uj 
24 ttHv 




(9) 


This has both a real and an imaginary part, logarithmi¬ 
cally divergent as function of A. Thus the DC conduc¬ 
tivity is zero, i.e. the material behaves like an insulator, 
qualitatively different from graphene that is a pseudo- 
dissipative metal. This already indicates that in 3D the 
Coulomb interaction is unscreened {2Qj. The dependence 
on the wave vector follows uniquely from the pseudo- rel¬ 
ativistic invariance of the free Weyl gas like in its 2D ana¬ 
logue, graphene |27]. Hence, using notations of Eq.(j4]), 
the electrodynamics is local with 

= uo) 


Interactions cause nonlocal electrodynamics 

Using the Coulomb interaction Hamiltonian within 
the tight binding model, Eq.(|6]), one obtains the correc¬ 
tions to first order in the effective interaction strength 
a = e 2 /kHv, where k denotes the dielectric constant of 
the background. While the corrections to either ar or 
aL are small, the relative correction to the difference 
a n i = aL — (Jt is actually dominant for k = 0 and can 
be considerable in the homogenous regime v 2 k 2 « uj 2 . 
While in free WSM a n i (uj, 0 ) was zero, now it becomes 
finite: 

(Jni (w, 0 ) = cr 0 (w) | ^ log + o (a 2 ) | . (11) 

Here uj is the frequency at which the renormalized value 
of the (renormalization group ”running”) Fermi veloc¬ 
ity v (uj) and the renormalized coupling a = a (uj) are 


2„,2 


A 2 v 


— v 2 k 2 


{Sij (uj 2 — v 2 k 2 ) + v 2 kikj } . 


( 8 ) 


r 


defined 20|, 22, 23]. Details of the calculation 


are given m 


Optics/Plasmonics within WSM 

Let us now develop the equations for the electromag¬ 
netic fields on the macroscopic scale in a WSM. We em¬ 
ploy the Landau-Lifshitz[30[ i3l| definition of the macro¬ 
scopic fields 


V D = 4 7T Pext ; VxE = -iB; (12) 

V • B = 0; V x B = + 47rJ ex ^j . 

The displacement field is D = E + 47rP, where the polar¬ 
ization vector is defined by P = J. The induced current 
Eq.([T]) is determined nonlocally by the conductivity ten¬ 
sor given in Eqs. m and ([IT]) . In particular, for the 
harmonic dependence E (r, t) = Eoe^ k ' r_a;t ) (and with¬ 
out external charges and currents) Gauss’ Law demands 
that either E is transversal, k • E =0, or else the disper¬ 
sion relation is plasmonic: 

47 r 

1 + i —<7 L (w,k) = 0. (13) 

UJ 

Ampere’s Law together with Faraday’s Law results in the 
condition that either E is longitudinal, k x E = 0, or the 
dispersion relation is transversal: 

. , c 2 k 2 

1 + i — a T (uj, k) = —— . (14) 

UJ UJ Z 

It is interesting to note that when uj « Av the con¬ 
ductivity, Eq. m is approximately purely imaginary and 
one obtains the dispersion relations for the longitudi¬ 
nal (” plasmon”, wavevector q) , and transversal (” light”, 
wavevector p), waves : 
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UJ 


1+1 + 


Na l A 2 v 2 \ 


-l 


24t r 


log 


7ra 0.6cj 2 

+ T ° s 


P 


2 



2+fiog(AV)_ \ 

1 + lo S (A 2 v 2 /w 2 ) J 


(15) 

(16) 


The results for q,p in units of uj/v as functions of uj/uj 
in the case of N = 4, v (uJ) = c/300, a (cJ) = 0.2, 
hvA = 5eV, (characteristic of Na^Bi and Cd^As^ for uj 
in the THz range) are presented in Fig. 3. One observes 
that the light wavevector p is real for all frequencies, 
that is the wave is nondissipating, see the left inset in 
Fig.3a. The plasmon wave vector q is real, according 
to Eq. ([T5]) . only above a certain threshold frequency uj p , 
see Fig.3b. One should emphasize that the ” plasma fre¬ 
quency” is that of a neutral plasma rather than that of 
the charged plasma in metals. However within the renor¬ 
malized perturbation theory used to derive the equation, 
the low frequency region might be unreliable. Above the 
threshold the two waves coexist and both wavevectors 


are real. 

In local materials equations Eq. (fl6|) and Eq. ([T5]) can¬ 
not be satisfied by the same wave vector: the difference 
of the equations is i^j-cr n i (cj, q m ) = —c 2 k 2 /uj 2 . Us¬ 
ing cr n /(u;,k) = — /? 0 (cj) /c 2 , see Eq. (flQ|) . this leads to 
A7riuj(3 0 (uj) = c 2 that generally cannot be satisfied. In 
WSM, on the other hand, cr n i (uj, 0) is finite and deter¬ 
mines the special value of the wave vector at which both 
transverse and longitudinal waves are the same. We fur¬ 
ther discuss this in the optical set up. Next we suggest 
two potential experimental setups that demonstrate the 
nonlocality of the WSM electrodynamics. 


Reflection and refraction on a dielectric - WSM 
interface 

Let’s consider the following setup, see Fig.4. An elec¬ 
tromagnetic wave with the wavevector k x > 0 ,k y = 0 
and k z < 0 is incident at the angle 0 with the normal on 
a semi-infinite WSM. The electric field with p - polariza¬ 
tion in vacuum is then described by a superposition of 
incoming and reflected waves: 


Er c = a 


fLo 

1 k ’ ’ 


k x 


k 


l J>(k x x+k z z) 




(■ k x x-k z z ) 


where A is the incoming and B the reflected amplitude, respectively. The field in the WSM is given by 


Ei = C T { —, 0 , ■ 

l P 


e i(k x x+p z z ) _|_ q 

I 


kx n Qz 1 

q q J 


,i(k x x-\-q z z ) 


(17) 


(18) 


where Ct and Cl are the transmitted transversal and 
longitudinal amplitudes. 

Employing the Landau-Lifshitz notation for the macro¬ 
scopic fields in Maxwell equations in dispersive media 
Ea. (fl2]) . used by Golubko v 13 1\ to resolve similar prob¬ 
lems as in exciton physics [32[, the list of the continuity 
conditions on the boundary is: 

(i) the normal component of D: 

rpvac Ani -rpwsm \ f 1 \ \ -rpwsm 

B z = —a zx E x + I 1 + — <r zz J A3 

(ii) the tangential component of E: E x ac = E x srn . 

(iii) the normal component of B: d x Ey ac = d x Ey Srn . 

(iv) the tangential component of B: d z E x ac —d x E^ ac = 

Q pjwsm _ Q pjwsm 

(v) In addition, since two different modes propagate in 
the bulk one needs a so-called ABC (additional boundary 
condition) [32, |33[. Assuming a ” sharp” interface (width 
smaller than the wavelength) the simplest ABC suffices: 


J z ( z = 0 ) = 0 ^ a zx E™ sm + a zz E™ srn = 0 . 

This reflects the fact that the electric current cannot es¬ 
cape the material. 

Applying these boundary conditions we obtain the fol¬ 
lowing reflection and transmission amplitudes: 

_ 1 - D _ 2k/p _ _ 2k^q (p 2 - fc 2 ) 

1 + D' T 1 + D' L kq z p 2 (\ + D) 
Ifpzqz - (p 2 - fc 2 ) kl 

k z q z p 2 ' 1 J 

When D = 1 the reflection coefficient vanishes and to¬ 
tal absorption occurs, thus turning the WSM opaque at 
frequency-dependent incident angles. 

The s - polarization does not generate a plasmon. The 
amplitudes are standard: r = (1 — p z /k z ) / (1 + p z /k z ) 
and t = 2/(1 + p z /k z ). In Fig.5a the amplitudes for the 
p-polarizations (solid) and the s-polarizations (dashed) 
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are presented. Fig.5b shows the vanishing of the p- 
polarization reflection coefficient at various incident an¬ 
gles in the THz range. 


Surface charging due to DC current 

Two possible experiments are suggested for low fre¬ 
quency transport in WSM. A Corbino geometry is em¬ 
ployed in order to avoid ambiguities related to bound¬ 
aries, see Figs.6 and 7. A slowly alternating flux is gen¬ 
erated by a thin solenoid inserted into the cylindrical 
aperture (radius R ) inside a sufficiently large WSM sam¬ 
ple. The fluxon’s magnetic field is linear in time (at times 
smaller than 1 /uS). Within that period of time the associ¬ 
ated electric fields in the aperture and the WSM are time 
independent. In the first setup the fluxon is concentrated 
at distance d from the center: <f> (r,t) = ^outS (r — dx). 
The solenoidal part of the electric field produced by the 
fluxon is 


E r l (r) = ^Je ij -b j , (20) 

where b = r — dx, i,j = x,y. It induces currents in the 
WSM that in turn generate charge distributions on the 
aperture surface, provided that the flux is not concentric 
with the cylinder axis. The surface charge density Q (0) 
is obtained from the condition that the normal compo¬ 
nent of the current on the surface on the semimetal side 
vanishes, n • J = 0. The irrotational part of the electric 
field is 


E" r (r) = L 



Q ( 0 ) 

(^2 + s 2 ) 3 / 2 ’ 


( 21 ) 


where the distance from a point on the surface to the 
observation point r is denoted by s = r—Rn'. The 
boundary condition, n • J = n- ( g t W° 1 + cri / E 2rr ) = 0, 
together with Gauss’s Law, n • E zrr =27rQ (0), leads to a 
dipole-like charge density. The details of the calculation 
are given in [34]. The surface charge density as a function 
of the polar angle 0 is 


QW 


— UJ&0 
CTL 


d sirup 

47r 2 (R 2 -2Rdcos(P + d 2 )' 


( 22 ) 


The dominant linear uj dependence of the conductivities 
cancels in the ratio ^ = 1 — 4^- log . The sur- 
face charge can be measured by gauging the electric field 
near the interior surface of the material. The deviation 
of cft /&l from unity (the value in the case of a local con¬ 
ductor) is an indicator of nonlocality of electrodynamics. 
Note that a central position of the fluxon would not in¬ 
duce the surface charges. The distribution of charge is 
given in Fig. 6b. 


for surface charging due to DC current, the calculation 
is more complex since induced charges appear not only 
on the aperture surface, but also on the 

For the second experimental setup, see Fig. 7 involve a 
Corbino disk with the fluxon in the center of a composite 
cylinder of WSM (0 < (p < n) and a usual ’’local” con¬ 
ductor (for example semiconductor) in the lower segment 
7T < (p < 27t; surface charges are now induced at the in¬ 
terfaces 0 = 0 and 7r between the two materials which in 
turn cause surface charges at the aperture. This results 
in an integral equation that fortunately can be solved 
exactly, see and results in 


QW 


a — (Jt <^4>o 1 + cos 0 

cf + cfl (2nf R S 1 - cos cf>' 


(23) 


The distribution of this surface charge density is shown 
in Fig.7. 


DISCUSSION 

The calculation of the leading Coulomb interaction 
effect on the electromagnetic response of 3D Weyl 
semimetal reveals that its macroscopic electrodynamics 
becomes nonlocal in a sense that the wave vector de¬ 
pendent AC conductivity tensor becomes nonanalytic at 
small wave vectors, see Eq. &■ The origin of nonlocality 
is a unique combination of the long - range (unscreened) 
Coulomb interactions like in dielectrics and the ultra - 
relativistic nature of the quasiparticles. The longitudi¬ 
nal conductivity gl (related by charge conservation to 
the electric susceptibility) and the transverse conductiv¬ 
ity gt are different in the long wave length limit and 
consequently the standard local Ohm’s law description 
does not apply. 

Several experimental signatures of the nonlocal contri¬ 
bution to conductivity, Eq. ([TT]) were pointed out. After 
reformulating the macroscopic electrodynamics in terms 
of the longitudinal and transverse fields with correspond¬ 
ing material parameters, the charging effect in DC trans¬ 
port in Corbino geometry was worked out. The set up 
avoids potential complications with conventional leads 
and contacts, while allowing for an unambiguous indica¬ 
tion of nonzero difference a n i = gl~gt- A more conven¬ 
tional set up is the measurement of optical response of 
the WSM that is also sensitive to the nonlocality. The p- 
polarized light at general frequency uo generates in these 
materials besides the transversal waves also the longitudi¬ 
nal one (that can be viewed as a charge density wave or a 
bulk plasmon). The propagation inside the WSM is only 
slightly attenuated, so that one can measure both the re¬ 
flection and two transmission amplitudes, all dependent 
on g n i . At a specific (material parameter dependent) 
frequency the two modes coincide, a phenomenon impos¬ 
sible in a local medium. The polarization dependence of 
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the light reflection (transmission) off a WSM slab is sen¬ 
sitive to both &L and Gp • Remarkably for any frequency 
there exists an incident angle where total absorption oc¬ 
curs, so that WSM is opaque. 

The applicability of perturbation theory in 3D Weyl 
semi-metals was quite recently addressed by several 
groups who used renormalization group and other non- 
perturbative methods like the random phase and large 
N approximations [Hi I35I . l36| . The latter shows pertur¬ 
bation theory to be reliable over a wide range of val¬ 
ues of the interaction strength a, up to a critical value 
a c « 14. In addition, typical values of the background 
dielectric constant in 3D WSM in the range k ~ 20 — 40 
have been measured [HI, which assures a below the criti¬ 
cal value of the interaction strength. In contrast, in the 
2D case, where a critical value a c = 0.78 was found [37[, 
the situation for perturbation theory is less fortunate 
since the measurement on samples substrated on boron 
nitride or on suspended samples provided values close to 

that m, mil. 
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• Fermi point 




b 


a 

FIG. 1: Spectrum of a charged metal contrasted to that of a Weyl semimetal. 

a. The Fermi level surface area is nonzero for a metal (or a band semimetal). The density of states at the Fermi level is finite 
(charged plasma) as is the density of charge carriers. 

b. In a Weyl semimetal at neutrality (Dirac) point the Fermi level crosses the spectrum at finite number of states, so that the 
Fermi level surface area is zero. The electron gas constitutes a neutral plasma. 




q 



FIG. 2: The tight binding model and Feynman diagrams for conductivity tensor. 

a. The cubic lattice. 

b. Propagators and Coulomb interaction. 

c. The leading order conductivity. 

d. The vertex correction. 

e. The self energy correction. 

f. The ”glasses” diagram. 


J 

j 
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FIG. 3: Two branches of excitations in neutral plasma of the isotropic Weyl semimetal. 

a. The transverse (red) and longitudinal (blue) wave vectors are monotonically increasing as function of frequency. At the 
intersection point at {qm,^m} the electromagnetic wave is monochromatic (lower right insert). 

b. The lower limiting (plasma) frequency at which the longitudinal wave appears as a function of the interaction strength a. 
The insert shows the upper limiting frequency. 



FIG. 4: Transmission and reflection of light in Weyl semimetal. 

a. The light beam is either reflected or transmitted inside the WSM, where it splits into a longitudinal (green) and a transverse 
(yellow) mode. 

b. The p-polarized light incident monochromatic beam with wave vector k incident at angle 0 is split in the WSM into 
longitudinal (green) and transverse (yellow) modes with wave vectors q and p, respectively. Directions of the electric fields are 
shown in blue. 




FIG. 5: Transmission and reflection of light in Weyl semimetal. 

a. Reflection and transmission amplitudes as function of the incident angle 9. The p-polarization amplitudes (solid): reflected 
(red), transmitted transversal (yellow) and longtudinal (green); the incident radiation is totally absorbed at the interface at 
angle 6 a • The s-polarization amplitudes are plotted dashed. 

b. The (p) reflection amplitude r as function of frequency for various incident angles displaying the total absorption. 
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FIG. 6: Transmission and reflection of light in Weyl semimetal. 

a. A thin solenoid inserted at distance d from the center of a hollow Weyl semimetal cylinder demonstrates the charging effect. 
The solenoid generates pulses of magnetic flux. During the pulse the internal surface of the aperture of radius R gets charged. 

b. Charging effect in a system combining WSM and a conductor. The upper segment is a WSM, while the lower is an ordinary 
’local’ conductor (a semiconductor or a metal). The solenoid is placed in the center of the aperture. Charges appear both on 
the internal surface and interfaces. The surface charge density [23] as function of polar angle with d/R = 1/2. 
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Part I 


Supplemental Information 1 

1 Calculation of the interaction corrections to 
conductivity of WSM 

1.1 The current density in the tight binding model 

The minimal substitution determines the coupling of the external electromag¬ 
netic field, described by vector potential A{ with electrons on the lattice: 

K mc [A (r, t)\= l ~Y, r n ,ic“Vfc£ +ai + he, (SI1(1)) 

n,i 

where the hopping integral T n ^ becomes 

r n ,j = 7exp dj(n + saj,t)| . (SI1(2)) 

Applying the charge symmetry transformation —)> e ZXn ^c^ the current den¬ 

sity on the link is 

J ni (t) = (SI1(3)) 

In linear response the current density operator is expanded up to the first order 
in A ni as J ni ( t ) = (t) + (t): 



- —c a t 

2 h n 

°f C i+ a, + h C 


2 

Pi 

Jti it) 

e va 

= 1 - 0 

2 ch 2 

/ Ai (n + sa, 
Js= 0 


Normalized operators in momentum space are = A/" -1 / 2 ^^e _2k n cg, where 

k 

the number of unit cells is AT = V/a 3 with V = a 3 1, where ki = sin (afcA 


There are eight Weyl points at which = 


= 0 inside the Brillouin zone. 


The calculation therefore was performed for TV = 8. 


1 






1.2 Matsubara Action and Feynman rules 

It is useful to represent the electron gas via the Matsubara action (r = it) in¬ 
volving the 3 + 1 Fermion field ^ and the static photon (auxiliary) field : 

( T ) (+ “ ie< + (0 - \ (d ( T ) n+a, ( T ) + ^ c ) 

+ 5 ( 2 ^n ( T ) - ^n+a; ( T ) “ ^n-a, (+ ( T ) 

(SI1(5» 

We are using convenient units where a = v = h = 1. The Feynman rules can be 
read from the Fourier transform 

^ = Yi \ V’L k-o-) (i - cos ( Pi ))<f> pu 

kco l PV ) p vi 

(SI1(6)) 

The fermion propagator is, 

G pu = (iw- p • a)" 1 = tU} ^ - , (SI1(7)) 

or + £p 

while the photon propagator is frequency independent, g puj = v p . The vertex is: 
ie^k-k'+p^w-u;'- These Feynman rules are used to calculate the linear response 
diagramatically. 

2 Direct calculation of the conductivity tensor 

2.1 General relation between conductivity tensor and the 
dielectric constant 

We calculate the transversal and longitudinal conductivities (Jt and <jl- From 
Eq.(3) we derive expression for two scalars formed from the conductivity tensor 



2<tt (<n, k) + ol (oj, k) = era (w) (SI1(8)) 

a L (oj, k) = lim (u>). 

Of course, due to charge conservation the transverse conductivity is proportional 
to the dielectric constant (Matsubara) [1] 

a L (oj,k) = +(l-e(w,fc)). (SI1(9» 
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2.2 Non-interacting WSM 

The ac conductivity of the non-interacting WSM on the lattice is obtained from 
the Kubo formula. The trace is (for positive cu) 

o n (hj) = ^ ^ ^ Tr \jj iGp V (J{ (G p v-\-u G p jy)} = 2e \ ^ t ^ ~ 

u w p W + 4 ^p) 

(SI1(10)) 

while the longitudinal conductivity is 

<ti(w) = lim-^^TV{k-aGp, I ,k-a(Gp )i , +u ,-G'p )! ,)}(SIl(ll)) 

p,zy 

2 2 \ CJ 1 

3 e ^ £ p(^ + 4 e 2) =3 aii - 

From here one deduces that = &L and thus no nonlocal conductivity arises 
in the non-interacting relativistic system, see general arguments in [27]. 

On the lattice one obtains 

4"> = 4"’ = ^ (- 5^ + 5§i log m) - (SI1(12)) 

This result is consistent with [20]. The dielectric constant in leading order is 

e (0) = e 2 (7 _ A log (w)^ . (SI1(13)) 

2.3 The interaction corrections: 

2.3.1 The trace in Eq.(SIl(8)) 

The three interaction corrections to the trace of the conductivity tensor are the 
self energy, the vertex and the ’glasses’ corrections, see Figs.2d-f, 


Ga, 




— E Tr K-p p • q G Pj „G q> pcr*G q jP+c^Gp^-j-^cr^] • 

^ pq up 

2p 4 

-E Tr [v q -pp • p Gp^G^pGp^aiGp^+vGi] ;(SI1 (14)) 

^ pq up 


UJ 


fk E TrGp^Gk+p. 






where the first term is the vertex, the second term the self energy and the third 
term the ’glasses’ diagram. The following shorthand is useful: 


k n 


sin ki ; ki = cos (ki ); 2 ki = sin (2 ki) ; Sk = 

n ^ n ^ n'} ~ ( ~ n ~ n ~ n 'j 

k\ , h > 2 , h j, k = , fc 2 , fc 3 j. 
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Calculating the trace and integrating over the internal frequencies v and p 
results in the following expressions for the three diagrams (symmetrized in p 
and q). The vertex 


e 4 N 


(jVer _ _ 


l ’q p C^i + ^2) 


id 288 qp £p£q (^2 _j_ 4 e 2 ) (^2 + 4^2 ) 


2 ’ 


(SI1(15)) 


with 


-1 = 


S 2 = 2 


the self energy, 


= -w 2 (4 (p • q) (p • q) - 2 p • 2 q| ; 

(P • q) (2p • 2q) + 12p • q - 16e 2 p q + 24q 3 • p 1 
h4(p-p) (q • q) ( p • q) — 8 (p • p) (q 3 • p) + p <—>• q j 


e 4 N _ v„_ 


= —E 

, . O/l 


q-p 


—3 ~r ^4 


id 24 qp £ p £ q y £ 2 (^2 _|_ 4 ^ 2 ^ 


(SI1(16)) 


=3 = lo 2 {-e 4 (p -q) + 2 e 2 (p 3 • q) + 2 e 2 (p • q) - (p 2 • p 2 ) (p • q)} ; 

S 4 = 4e 2 p 4 (p • q) + 2e 2 (p 3 • q) - 2e 2 (p • q) - 3 (p 2 • p 2 ) (p • q)} , 

and the ” glasses” diagram that proportional to the square of the number of the 
Dirac points points, 


g i _ e 4 4t tN p 2 y Px \ 

<*> 3 ytel^ + Ael)) 


(SI1(17)) 


The integration over the momenta is performed with the sum of the vertex 
and the self energy diagrams, which ensures convergence. The integral for the 
glasses diagram vanishes. 


41 = -e 4w Q + log (w) + T- i 0 g 2 ( w )^ (SI1 (18)) 

which expressed in physical units and physical frequency is 


(i) _ 


a • • = 


e 2 aN 

—-- Ud 

hvS 


1.308 

2 




1.308 , 
-log 

7 r 




2.3.2 The longitudinal conductivity via the dielectric constant, Eq.(SIl(9)): 

The general expressions for the dielectric function via p — p correlator can be 
written as 5 = 1 + e 2 £ x + e 4 £ 2 with 

e 1 = - e 2 v k Tr G P)I/ G p+ k,„+u> (SI1(18)) 

pu 


4 



and s 2 = e ver + £ se + e gl where 


g .ver 


£ 


se 


£ al 


& 'Y2 TrG^p-|-k,zy+u;^p,^^q^p+q,p^p+q+k,p+u; 

pq vp 

2e v k ^q—pTrG^p ?Z yG^q j pG f p ?Z yG^p-|-k,a;+^ 

pqzyp 


^ k Tr£G p ,pG p+k 


, 0 ;+^ 


(SI1(19)) 


The terms are expanded in powers of the trace is calculated and integrated 
over the internal frequencies. All the k° and k 1 terms vanish after integration 
over the internal frequencies. 

The relevant terms (symmetrized in p and q): 


e 4 7 T iV Vq- p 

a 




( 3 ~ £ p)(^ 2 ~ 4£ p)+ 2 (^ 2 + l 2g p) _ [3u; 2 +2Qgp](2p) 

2( W 2 +4 £ |) 2 e | 

_ (^ 2 + 12£p)p 3 -q 

(o; 2 +4£ 3 ) 2 £ : 


8ep(w 2 +4£|) 2 * ^ ^ 


+ P 


(SI1(20)) 


e 4 TrN 

48 


£■ 


£p£q ( w2 + 4 £p) ( w2 + 4 ^q) 

/ [ 4 e p £ q - w 2 (p • q)] ( 2 p • 2 q) - 2 e 2 { 2 e 2 ( 2 p • 2 q) - w 2 2 q • (q * p 
1 +4e 2e q [4 (P • q) - w 2 ] (p • q) + p 4—» q 


(SI1(21)) 


s° l = 


eW 2 j ^ ~£ 2 + |(P 2 ■ P 2 ) + j|jp 
£p (w 2 + 4e2) 


(SI1(22)) 


The numerical integration over the internal momenta results in (expressed 
in physical units, Matsubara frequency) 


_ f N [-0-055 + 0.023 log (if) log 2 (^)] 

6 ^ a l +(I) 2 [i-bog (^)] 2 

and in physical frequency: 


(SI1(23)) 



N 


-0.055 + 0.022 log (^) - log 2 (^) + ^ + * (0.036 - ± log (^)) 


+&f 


(SH(24)) 
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3 Supplemental Information 2. Surface charg¬ 
ing 

3.1 Corbino geometry with a fluxon out of the center 

Assume that the disk is thick enough so that the system is z translational 
invariant. A thin fluxon placed at distance d on the x axis. The flux therefore 
can be represented by (r) = atS (b); b = r — dx; see Fig.6. This produces the 
rotational electric field 


E ! 0l (r) = ^e ij b r (SI2(1)) 

Surface charges are induced at the aperture with density Q (0). The density 
is obtained from the condition that the normal component of the current on the 
aperture surface on the WSM side vanishes, n • J = 0. The current is caused 
by the solenoidal field with conductivity gt and the irrotational field (induced 
by the surface charges) with conductivity gl • The irrational field (denoting the 
distance from a point at the surface to the observation point r by s = r— Rn') 
is: 


E? 




(‘n-Rrii ) Q (4>) 


(SI2(2)) 


Inside the semi-metal near the surface one has for its normal component (see 
Appendix) 


n • E irr =27T<2 (4>). (SI2(3)) 

The above condition is then 
oltl ■ 

cI t --- 2 £i i (^ n - (0) = 0, (SI2(4)) 

27r|i2n — dk\ 

so that one obtains Eq.(22) 

3.2 The WSM - (Ohmic) conductor circuit 

The fiuxon (in the -\-z direction) is located in the center of a composite cylinder 
of WSM with conductivities gt,<Jl in the ’upper’ (u) half (0 < < n) and a 

local conductor with Ohmic conductivity g in the lower (d) segment i t < (j) < 2i r. 
Surface charges are induced at the aperture surface and at in the interface 
between the two materials. Their densities are determined by the condition 
that the normal component of the current on the aperture surface vanishes. 
Only the irrotational field is involved which has contributions from the two 
interfaces and the surface charges. The contribution from the circular surface 
is dominated by the local charge (see Appendix). E (x) denotes the interface 
charge density and Q ((/>) the aperture surface charge density. 
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In the upper segment, the WSM, the condition at the boundary is then (in 
the lower segment the condition is trivially fulfilled) 

r (r — {x, 0}) S (x) 

R (r-{a;,0}) 2 
(SI2(5)) 

where the 2 : integration was already performed, as in Eq.(SI 2 ( 2 )). Changing 
x —x in the second integral, one gets 

2 f°° / R 2 — xRcos(t) R 2 + xR cos </> \ 

R J x= r X \ R 2 — 2 Rx cos (j) + x 2 R 2 + 2 Rx cos <p + x 2 ) ^ 

(S12(6)) 

Consider the right interface. Across the interface the current is continuous. 

Both above and below the interface the current consists of the solenoidal (trans¬ 
verse) and two irrotational (longitudinal) contributions from the interface and 
the aperture surface. The upper is 


J r =a L El rr (R^) =a L 


rOO 

' 

J x=R 


r (r-{x,0})S(x) 
R (r-{a;,0}) 2 


r-R 


J * (l ' 0) = “Sf + (a]) + 2i7l L„ (SI2(7)) 


while the lower is 


4 (*> 0) = ~T~ - (j27rS (*) + 2(7 r 2 V ( SI2 ( 8 )) 

^ 2tyx Jj ) f = 0 x 2 + R z — 2xRcos<p 


Equating the two, one obtains: 


ex f 2n 

(a - <j t ) - -h 27T (cr + cr L ) £ (x) - 2 (a - a L ) / 

Z 7 TX J (k' = 0 


R sin (f)'Q (0') 


= 0 


/ 0/ = 0 x 2 + R 2 — 2xRcos(j)' 

(SI2(9)) 

These two integral equations can be solved analytically. Rewriting Eq.(SI2(8)) 


as 


a f 27T 

kt~ -b 2i rE (x) = 2rlR / 

2tt x C =0 


sin cj)Q (0) 


=0 x 2 + R 2 — 2xR cos 0 ’ 


(SI2(10)) 


with and kr = and substituting Eq.(SI3(7)) into the latter 

gives 


—+27rS M ^ C sin ^ /°° s ^ f os0 

27tx 7 r 0 x 2 + R 2 — 2xR cos 4> J x '=r \R 2 + x' 2 — 2Rx' cos </> 

(SI 2 ( 11 )) 

The integration over </> on the right hand side vanishes leading to 

(SI3(12)) 

(27r) x 

This combined with Eq.(SI2(7)) leads to the result Eq.(23). 


2ttQ (</>) | r=jR = 0, 


R 2 + x'R cos 
R 2 + x /2 + 2Rx' cos 0 
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4 Appendix. The electric field at a circular aper¬ 
ture with an overall neutral surface charge. 

A conducting cylinder with a circular aperture (radius R ) is charged locally with 
surface charge density Q (0), albeit overall neutral: Q (</>) = 0. We are in¬ 

terested in the electric field just inside the conductor at R-j- A. The vector from a 
point on the aperture surface to the observation point r = (R + A) {cos y, sin y} 
is s = r— Rn', n' = {cos </>, sin (/)}. The (irrotational) field is given in SI2(2). 

The normal to the aperture component just inside the metal is 


—-E zrr = 2 
R 


r2ir 

Jcb =o r2 


r 2 -Rn'r 
+ R 2 — 2Ln'r 


Q (</>) = 2 


r2n 

J(b =0 


After some algebraic manipulations one gets 


(L + A) 2 —L (L + A) cos (<t>-x) 

=o L 2 + (L + A) 2 -2 L (L + A) cos (</> - X ) 
(SI2(13)) 




—-E* rr = 2 Q 
R 


(x) f 

J6 


(L + A) z —L (L + A) cos (4>) o f^~ e (L + A) 2 —L (L + A) cos 4> 

-+ 2 / r o ,„,,' 2 rT 77 — r ~ - : Q ( 4 >- 

J 6= 


,=-. e L 2 + (L + A) —2L (L + A) cos (0) L 2 + (L + A) -2L (L + A) cos 


The integral in the first term is 


/ 


(£ + A) 2 -i (L + A) cos (0) 
=- e £ 2 + (£ + A) 2 —2L (L + A) cos (0) 


= 7T + 6 — 2 arctan 


(SI2(14)) 
A cot (e/2) 


2L + A ’ 
(SI2(15)) 

and in the A -A 0 limit i t + e. In this limit the second integral is simply 

f£‘Q(4 + x). 

Taking now e -A 0 and using neutrality one has the result that the contri¬ 
bution from the circular surface is dominated by the local charge: 


-■E irr =2nQ( X ). 


(SI2(16)) 
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